Bounds on neutrino magnetic moment tensor from solar neutrinos 



Anjan S. Joshipura and Subhendra Mohanty 
Physical Research Laboratory, Navrangpura, Ahmedabad - 380 009, India 

Solar neutrinos with non-zero magnetic moments will contribute to the electron scattering rates 
in the Super-Kamiokande experiment. The magnetic moment scattering events in Super-K can 
be accommodated in the standard VO or MSW solutions by a change of the parameter space of 
mass square difference and mixing angle-but the shifted neutrino parameters obtained from Super-K 
will (for some values of neutrino magnetic moments) become incompatible with the fits from SNO, 
Gallium and Chlorine experiments. We compute the upper bounds on the Dirac and Majorana 
magnetic moments of solar neutrinos by simultaneously fitting all the observed solar neutrino rates. 
The bounds the magnetic moment matrix elements are of the order of 10" 10 /is. 
1. Introduction 

Solar neutrino experiments with Chlorine Gallium [||, and water cerenkov detectors ||, show that there 
is a deficit of v e at the Earth compared to the predictions ||] of the standard solar model (SSM). It is commonly 
accepted that vacuum oscillations (VO) or matter induced MSW conversions, with large mixing angle (LMA) or small 
mixing angle (SMA), can account for the deficit of solar neutrinos. The available experimental results allow not only 
^ ■ a test of the neutrino oscillation hypothesis but also offer possibilities of constraining new physics, e.g. neutrino 
magnetic moment [^|,0j9), neutrino decay |Io[ , flavour changing neutral currents, etc. 

Neutrinos with non-zero magnetic moments will contribute to the ve~ scattering rates at Super-K. The electron 
scattering rate at Super-K (defined as the ratio of average number of electron scatterings observed and the theoretical 
■ prediction of v e e~ elastic scattering in SSM) is given by 

m : 1 

O ■ R SK = — ((a e P ee ) + (a M , T (l - P ee )) + (a mag )) (1) 

-xt- : 

Here the first two terms denote the standard v e e and v^. T e scattering processes respectively , P ee being the v e 
survival probability. The last term in eq.(Q) corresponds to the scattering events — > Vj&~ due to non-zero 
neutrino magnetic moments. Unlike, Super-K and SNO, the other solar detectors do not probe the neutral current 
interactions and are thus not sensitive to the magnetic moment scattering. The previous bounds on neutrino magnetic 
moment from solar neutrinos have made use of the only Super-K results ||[7| or have combined Super-K and SNO 
O ,' results Q together. The other experiments however indirectly influence the bound on magnetic moment through a 
combined chi-square analysis. We use results from all experiments here concentrating mainly on rates to put bounds 
' on elements of the magnetic moment tensor fiij defined in the neutrino mass basis. 
^ , The magnetic moment scatterings are helicity flipping processes and are more appropriately described in the mass 
basis unlike the non-flip weak interaction processes which are described in the flavor basis. The electron neutrinos 
produced at the Sun arrive as a linear combination of mass eigenstates at the detector. The probability amplitude 
of the Vi mass eigenstate at the detector is denoted by Ai{L) and i^e - — > Vje~ scattering amplitude is proportional 
to fiij . It is assumed that prior to scattering the neutrinos are a linear combination of different mass eigenstates, 
but after scattering the neutrino beam is a incoherent sum of different mass states. A field theoretic justification of 
this is given in j^] . The probability amplitude of the Vi mass eigenstate at the detector is denoted by Ai (L) and the 
z^e~ — * v j e ~ scattering amplitude is proportional to /Ay. The magnetic moment differential scattering cross section 
at the detector is then given by 

dcJmag 2 A ^ \ fn \ 

T is the recoil energy of the electrons, E v the initial neutrino energy and the magnetic moments fiij are in units of 
Bohr magnetron. The effective magnetic moment fi e ff ca n be written in terms of the intrinsic neutrino magnetic 
moments /ly as 

j i 

The magnetic moment tensor /xy is antisymmetric if neutrinos are Majorana particles. For Dirac neutrinos, fiij are 
arbitrary complex numbers but in special case of CP conserving theory, they are real and symmetric f|. We shall 
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assume this to be the case. The amplitude Ail of the z/j mass eigenstate at the detector depends upon the oscillation 
model which solves the solar neutrino problem. 

We consider the case of vacuum oscillations, and matter induced MSW conversions in turn and derive bounds on 
Hij for each of these cases separately. 



2. Vacuum Oscillations 



Neutrinos are produced in the solar core as the v e flavor eigenstate, and the amplitude for detecting them at a 
distance L as a mass eigenstate is given by 









H 








ML) = U el e^ L . (4) 

The elements of mixing matrix U are constrained by atmospheric neutrino and CHOOZ experiments to be of the 
following form: 

I v e \ I ci si \ / v\ \ 

(5) 

We have approximated U e 3 ~ in order to account for the negative results obtained in v e disappearance experiment 
at CHOOZ jll). The small value of t/ e 3 (taken here as zero) also prevents the mixing of the third neutrino in vacuum. 
In addition to the above mixing one also needs 8m\ 2 = m\ — to 2 « 10~ 4 — 10~ n eV 2 and 5m\^ = to 2 — to 2 , sa 10~ 3 
eV 2 to account for the solar and atmospheric Jl2| neutrino problems respectively. The required value of 8m\^ is 
much larger than the value of the effective mass square 2^2GpEN e of the electron neutrino at the solar core. This 
suppresses mixing of the third neutrino in matter as well |l5} ] which decouples from the rest in case of the MSW 
solution to the solar neutrino problem. 

Using the mixing matrix (j^) the v e survival probability in vacuum from Sun to Earth is given by the expressions 

P ee = c{ + s\ +2(c lSl ) 2 cos(^|^) (6) 



500.3 s is the Earth-Sun distance. 



Using the mixing matrix elements from eq. (o) , we find that the expression for /i 2 f f for VO is given by 



3 V i 



3 

= c i (Mn + M? 2 + Mis) + 4 (f4i + M22 + M23) 

+ 2cisi (/X11/X21 + M12M22 + cos(— — — ) (7) 

lhi v 

This is the three generation generalization of the expression for (/i 2 j^)yo first derived by Beacom and Vogel The 
1/3 state can be produced after the electron scattering by v\ and V2 in the solar neutrino beam, therefore the transition 
moments [in and /i23 appear in the expression above. However since the solar neutrino beam does not contain v$, 
there is no [i^i dependence in fi e ff and /133 cannot be constrained by solar neutrino experiments. 

Vacuum solution requires 5m\ 2 ~ 10~ 10 eV 2 . The oscillatory term in cq.(Q) does not get averaged out in this case. 
This term by itself is not positive definite and prevents obtaining bounds on all the moments appearing in eq. (Q) . We 
thus need to make simplifying assumptions. For the case of Dirac neutrinos we assume (as in the Particle Data Book) 
that the dominant non-zero elements are diagonal fiij = SijfXi. With this assumption, eq.(Q) simplifies to the form 

2 2 2,22 /a\ 

Veff = C l/»1 + SlM2- ( 8 ) 

Since each term ^ is positive definite we can put upper bounds on n\ and by taking them one at a time. 

For the case of Majorana neutrinos the magnetic moment matrix fj,ij is anti-symmetric. We set /in = ^22 = and 
find that /i e f / can be given an upper and lower bounds 

M? 2 + (ci|/*ia| " si|m 23 |) 2 < Me// < M?2 + OilMial + si|m 23 |) 2 - (9) 
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From this we see that A*^// can be written as a sum of two positive definite quantities and therefore one can take each 
term separately in order to put an upper bound on (j,± 2 - To put bounds on /J43 and ^23 we need to further assume 
that there is no cancellation between the two terms in the the expression (ci I/X13 1 — Si |yu,23 1) 2 i.e either I/X23I or | >Lt 13 1 
is much larger than the other. 

3. MSW conversion 

If the flavor conversion of neutrinos in the Sun is by MSW mechanism then we have different expressions for fj, e ff. 
As discussed already, only two energy eigenstates v\ and v-i corresponding to the lighter neutrinos participate in the 
MSW conversion. In the core of the Sun E\ > E2, there is a level crossing at the resonance point after which E% > E\ 
i.e. mi > mi in vacuum. The probability of v e — > v\ just after level crossing is 

Pi = Pjs 2 m + (1 - Pj)c 2 m (10) 

Here, the first term stands for v e going to V2 by mixing with probability and then jumping to v\ with probability 
Pj at the level crossing. The second term means that v e goes to v\ by mixing and then does not jump to with 
probability (1 — Pj) at the level crossing. The probability of v e — > just after level crossing is 

P 2 = (1 - Pi) = (1 - Pj)s 2 m + Pjc 2 m . (11) 

In the formulas above, the Landau- Zener jump probability is given by 

_ exp(-bsj /E v ) - cxp(-6/g„) 
3 ~ 1 - exp(-b/E v ) [iZ) 

n ( ^-10»f^MeV (13) 



4 \\A/A\ re J VeV 

and A = 2^/2E v GpN e . The mixing angle in matter in the Sun is given in terms of the vacuum mixing angle by the 
expression 

cosZUm - (1-27,(1- 2 S \) + v 2 y/ 2 1 ' 

n= S ™L = 6.6X10-^ b_ (15) 

We can now use the matrix eq.(|5|) to find the survival probability of v e at Earth, 

P ee = c\P x + s 2 P 2 (16) 

Unlike in the vacuum case, 2 ^% is much smaller than the Sun-Earth distance in case of the MSW conversion. Thus 
phases acquired in oscillations average out juj. We can thus drop the interference terms of the type A\{L)A2{L)* in 
the expression for /i e // m ec f(||)- Substituting for |^l 1 (L)| 2 = Pi and |^l 2 (i')| 2 = P2 an d we then get 



(^eff)MSW = Y. P ^l 



= Pi (Mil + M12 + M13) + p 2 (f4i + M22 + M23) 

= rt 2 + Pi( Mll + 14 3 ) + (1 - Pi)04 + M 2 a) (17) 

The expression for (/i 2 ^ j)msw was first derived for the two generations case in Q. Since, (/i^ j)msw contains a sum 
of positive definite terms we can put upper bounds on each term taken one at a time, using the expression eq.([Tc|) for 
Pi. In this case, we are able to constrain all fiij except ^33. 

4. Experimental rates and bounds on magnetic moments 

The electron scattering reaction in Super-K can be written as 
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/ dE v <T Vc $ B Pee + J dE v a v ^ B (l - Pee) + J dE (^) $ B 

J dE u a Ue <I>B 

The magnetic moment scattering cross section a mag (E u ) is obtained by integrating (w.r.t the electron recoil energy 
T )the differential cross section || folded with a detector response function taken from p5| . The 8 B v e flux is taken 
from H and is given below. Only the Super-K electron elastic scattering rate depends on the neutrino magnetic 
moments. 

The charge-current deuterium dissociation reaction rate at SNO normalized to its SSM value is given by 

cc jdE u a GC ^Pee , . 

RSN ° ~ fdE„* cc t> B (19) 

The v e e~ and v^. T e~ elastic scattering cross section after folding with the detector response function are tabulated 
in ref . Jl6| . The deuterium dissociation cross section is taken from JlTj . 

The rates of neutrino capture in the Chlorine and Gallium experiments can be written as 



E 



/ dE v §i<j a P e , 



R a = !=ggg (20) 

2^ J dE„<$>i(T a 

i— pp. Be, B 

where the subscript a — Ga, CI denotes the experiment and i — pp, 7 Be, 8 B denotes the type of neutrino flux from 
the Sun. The spectra of the pp and 8 B neutrinos can be fitted with the analytical functions, 



<S> pp = (5.95 x 10 10 )[193.9(0.931 - £„)((0.931 - E v f - 0.261) 1/2 £ 2 ] 
$ B = (5.05 x 10 6 )[8.52 x 10~ 6 (15.1 - E v ) x75 El] 

$ Bo = (4.77 x 10 9 )[5{E V - 0.862)] (21) 

where the neutrino fluxes are in units of cm _2 s _1 and E v is in MeV. The first brackets in ( |2l] ) give the total flux of 
neutrinos from the pp,Be, and B reactions and are taken from BP2000 ||, and the square brackets give the spectral 
shape @. 

The Ga experiments can detect all three types of neutrino fluxes and the neutrino absorbtion cross section of Ga 
is given in |^9|. The Chlorine experiment threshold is higher (0.8 MeV) and it detects only the Be and B neutrinos; 
we take the absorption cross section with the tables from ref. pjj . 

Using the flux spectrum in equation ( |2l| ) and the cross sections |T^,^,|l6|,^7) we can calculate the theoretical rates 
for MSW or VO conversion probabilities as a function of the three unknown parameters: the magnetic moment 
matrix elements /Lt^, 5m\ 2 and the vacuum mixing angle Q\. The experimental rates i?° xpt with one-sigma combined 
(statistical and systematic) experimental errors A Q are as follows |^-^] : 

R c ^ = 0.335 ± 0.029 
R™^ = 0.584 ± 0.039 
i?^ pt = 0.459 ± 0.017 

#SNO-CC = °- 347 ± °- 029 ( 22 ) 

From the theoretical R a , and the experimental i?° xpt , we compute the total x 2 for all experiments, defined as 

/ r? Dcxpt^2 

X 2 - E — ~^ ] ■ (23) 

a=Cl,Ga,SK,SNO a 

Setting the parameter fi e ff = we reproduce the standard contours of VO (Fig. 1.) and MSW (Fig. 2) solutions 
shown as continues curves. For ls e ff = the global minimum of \ 2 occurs in the vacuum region with Xmin ~ O- 3 - 
We then take a non-zero /z e / / and plot the regions (shown as dashed curves) in 5m\ 2 — tan 2 plane corresponding 
to Xmin + 4-6 (90% CL) both for VO (Fig.l) and MSW ( Fig. 2) solutions. These figures correspond to a constant 
value for f. It is found that the allowed regions shrinks as we increase the value of /i e f / and it disappears for some 
specific /jb e ff which is taken to be the 90% CL bound on /Lt e //. We consider different cases following this procedure. 
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In case of VO and Dirac neutrino, we choose m 2 // equal to n\c\ and i3 L s\ (see eq.(||)) respectively. The bounds on 
Mi. 2 obtained using the above procedure are given by (at 90%C.L.) 

Mi < 4.5 x 1CT 1c Vb 

H2 < 7.1 x 10- 10 mb (24) 

In case of the Majorana neutrinos we use eq.(Q) and substitute Me// by M12 i c i IM13 1 2 or -s 1 1 M23 1 2 leading to the 90% 
CL bounds 

IM12I <3.7x 10- 10 mb , 

|Mis| <4.5x 10- 10 mb (*/ |M2 3 |«|Mi3|) 

| M23 | < 7.1 x 10- 1( Vb (if |/ii3|«|M23|) (25) 

Next we take up the MSW solution. We find that the LMA and SMA allowed parameter space disappears for 
Me// = 3.7 which translates to the bound on Mi2- Next we take Me// = (1 — -Pl)m 2 an( i nn( i that the parameter space 
vanishes when Me// = 4.3. This translates to the bounds on I/X22 1 , M23 1 ■ Finally we take Me// = ^lM 2 an d find that 
the SMA region disappears for Me// = 7.8 and the LMA region disappears at Me// = 20.5 which translates into the 
bounds on |mii|> Misl- These bounds are summarized below: 

IM12I < 3.8 x 10- 10 mb 
IM22UM23I <4.3x 10- 10 mb 

|mu|, |Mis| < 7.8 x 10" 10 mb (SMA) ; 20.5 x 10 _1 °mb (LMA)) (26) 

In ( |26j ) , the bounds on the diagonal elements of My ar e relevant for Dirac neutrinos only and the bounds on the 
off-diagonal elements apply to both Dirac and Majorana neutrinos. The bound on |mii|, ImisI obtained in case of the 
LMA solution is weakest of all bounds. This is due to the fact that these quantities are multiplied by Pi in eq.(|l7|) 
which is very small for adiabatic transition occurring in the LMA region. 

While the neutrino magnetic moments affect the rates of only the Super-K experiment, a combined analysis of all 
experiments was essential in establishing the bounds given in eqs.( p4| - |26] ). This can be seen explicitly in Fig. 3. and 
Fig. 4 for the VO case and in Fig. 5. and Fig. 6. for the MSW case. In Fig. 3. we plot the 1.64a allowed regions 
for Rsk (continuous line) and Rsno (dotted line) by setting Me// = 0. We see a large overlap between the SK and 
SNO allowed regions for Me// = 0. Next we repeat the same plots by taking Me// = 3.8. We see that the SK allowed 
parameter shifts so that there is no more any overlap with the SNO allowed region. Therefore a combination of SNO 
and SK result is crucial in disallowing Me// > 3.8. 

Similarly in Fig. 5. and Fig. 6. we plot the 1.64a allowed regions for Me// = and Me// = 4.0 respectively. In 
Fig. 5. we see overlap between SK allowed regions (between continuous curves), SNO (between dashed curves) and 
Gallium (between dotted curves) in the SMA and the LMA regions. In Fig. 6. we see that when Me// = 4.0 the SK 
allowed parameter space shifts so that there is no common allowed region for the three experiments thereby ruling 
out values of Me// > 4.0. 

In summary, we have used solar neutrino rates to constrain possible magnetic moment coupling of neutrinos to 
photons. The bounds depend upon specific solution of the solar neutrino problem. If MSW mechanism is mainly 
responsible for the solar neutrino conversion then it is possible to constrain all elements of magnetic moment operator 
Mij except for M33- For the VO solution one can obtain bounds (eq.(|2j)) on two of the diagonal moments assuming a 
diagonal Mij in case of the Dirac neutrinos. For Majorana neutrinos, we obtain bounds on M12 > M13 an d M23 given in 
eq-dH). 
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FIG. 1. The 90% allowed parameter space for VO , with fi e ff = (enclosed by continuous curve) and for fj, e ff = 3 (enclosed 
by dashed curve) . 
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FIG. 2. The 90% allowed parameter space for MSW conversion , with fj, e ff = (enclosed by continuous curve) and for 
fi e ff = 3 (enclosed by dashed curve). 
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FIG. 3. The 1.64a VO allowed regions for Rsno ( enclosed by dashed curves) and Rsk ( enclosed by continuous curves) for 

fleff = 0. 
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FIG. 4. The 1.64cr VO allowed regions for Rsno ( enclosed by dashed curves) and Rsk ( enclosed by continuous curves) for 

fleff = 3.8. 
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FIG. 5. The 1.64er MSW allowed regions for Rsno ( enclosed by dashed curves) and Rsk ( enclosed by continuous curves) 
and Rgo, (enclosed between dotted curves) for p e // = 0. There is ovelap in the SMA and LMA regions. 




FIG. 6. The 1.64a MSW allowed regions for Rsno ( enclosed by dashed curves) and Rsk ( enclosed by continuous curves) 
and Roa (enclosed between dotted curves) for (j, e ff — 4. The SK allowed region shifts such that there is no ovelap between the 
three. 
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